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Abstract 


In  this  report  we  study  the  effects  of  the  non-linearities 
in  the  dispersion  and  attenuation  of  longitudinal  waves  in  a  plasma. 
We  first  work  with  the  hydrodynamic  equations  and  show  that  to  a 
first  approximation  non-linearities  affect  the  dispersion  equation. 
We  then  work  with  the  Boltzmann- Vlasov  equation  and  indicate  a 
general  method  for  the  study  of  the  effects  of  non-linearities.  Be¬ 
cause  of  mathematical  difficulties,  we  develop  a  hybrid  formulation 
which  is  based  on  both  the  hydrodynamic  and  the  Boltzmann-Vlasov 
equations  and  indicate  that  to  a  first  approximation  non-linearities 
affect  both  the  dispersion  equation  and  the  attenuation  in  a  plasma. 
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I. 


Introduction 


By  linearizing  the  Boltzmann- Vlasov  (B-V)  or  the  hydrodynamic  (H) 
equations  that  describe  the  longitudinal  properties  of  a  plasma,  in  effect  we 
disregard  interactions  between  modes  of  vibration  in  the  system.  If  accounted 
for,  these  interactions  could  lead  to  attenuation  of  an  excitation  (density  wave) 
in  the  plasma  in  addition  to  the  attenuation  due  to  collisions  or  to  Landau  damp¬ 
ing;  they  could  further  lead  to  a  modification  of  the  plasma  dispersion  relation 
which,  in  the  linearized  case,  is  given  by 

=  cOp  +  (3RT/m)  ;  (B-V)  (1) 

w  =  ojp  ;  (H)  (2) 

The  added  attenuation  can  be  attributed  to  the  energy  lost  by  a  propagating  den¬ 
sity  wave  for  the  excitation  of  harmonics.  On  the  other  hand,  the  modification 
of  the  dispersion  relation  can  be  attributed  to  a  correction  of  the  longitudinal 
dielectric  coefficient  of  the  plasma  because  of  non-linearities. 

Some  attempts  have  been  made  in  the  last  few  years  to  study  the  non¬ 
linear  effects  in  a  plasma  in  equilibrium  described  by  the  Boltzmann- Vlasov 
(1  2  3) 

equation  ’  ’  .  The  little  progress  made  consists  perhaps  in  the  realization 

of  the  mathematical  and  conceptual  difficulties,  but  has  added  no  new  insight 
to  the  problem.  In  the  present  paper  we  make  a  renewed  effort  to  understand 
the  problem  by  attempting  to  describe  the  non-linear  effects  as  a  scattering 
process.  In  Sec.  II  we  limit  the  discussion  to  the  hydrodynamic  model  in  order 
to  keep  a  simple  mathematical  structure.  We  are  aware  that  the  hydrodynamic 
model  is  inadequate  for  large  wave  numbers  k  because  it  does  not  account  for 


the  interaction  of  the  collective  mode  with  the  free  particles  (leading  to  Landau 
damping).  Although  this  shortcoming  does  not  affect  the  validity  of  our  method, 
it  does  cast  some  skepticism  on  the  accuracy  of  the  results.  In  principle,  the 
method  can  be  generalized  to  treat  the  Vlasov  plasma.  This  generalization  is 
developed  in  Sec.  III. 
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II. _ Non-Linear  Hydrodynamic  Equations  for  a  Plasma 

We  use  the  following  set  of  equations  to  study  the  non-linearities  in 

(4) 

a  homogeneous  model  of  a  plasma  ; 


(a)  Continuity  Equation; 


+  V-  •  [n{xt)  V  (xt)J  =  0  , 


X 


(3) 


(b)  Equation  of  Motion; 


8V(xt) 


+  [V(xt)  •  V-]  V  (xt)  +  ^  V<j>(xt)  =  0  , 

01  X  m 


(4) 


(c)  Poisson  Equation; 


V^<f(xt)  =  -47reSn  (xt) 


(5) 


The  charge  of  the  electron  has  been  taken  as  -e.  and  the  mean  electron  density 
n^  is  assumed  to  be  neutralized  by  a  positive  sea  of  ions,  ^n(xt)  is  the  excess 
density  at  a  point  x  and  time  t  and  therefore 


n(xt)  =  n  +  Sn(xt)  =  n  {1  +  J  =  n  {1  +  i1j(xi 


o'  n  '  o 
o 


t)} 


In  terms  of  i}j(xt),  Eqs.  (3)  and  (5)  become 


3 


+  V-  •  V(xt)  +  V-  •  [ijj(xt)  V  (xt)]  =  0 

A  ^ 


(3a) 


9iJj(xt) 

at 

2  —  — 
V  <j)(xt)  =  -47r  en^ip(xt) 


(5a) 


By  taking  the  Fourier  transforms  in  time  and  space  of  Eqs.  (3a),  (4),  and  (5a), 
we  obtain 


-icj  4;(k<j)  +  ik  *  V(kci))  +  ik  •  J  V  (k'w  ')4i(lc-E',  co  -co  ')  =  0  (6) 

(27r)^ 


-iwV(kw)  +  ik—  (j>  (ko))  +  i/(k-k')  •  V  (Pw ')  V  (E-k',  w  -  w')  =  0  ( 

{ZTrr 

-k^<j)  (kw)  +  47r  en^ijj(k  w)  =  0  ,  (8) 

where  the  quantities  in  the  integrals  represent  the  non-linear  contributions. 

The  linearized  version  of  the  hydrodynamic  equations  have  the  form 


-iwi|^(kw)  +  ik  •  V(kco)  =  0 

(9) 

-iwV(kw)  +  ik ^  (kw)  =  0 

(10) 

-k^<j)(Ew)  +  4  TT  en^i|j(^‘^)  ~  0 

(11) 
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and  yield  after  some  trivial  substitutions 


^2 

k  {1 - 4)  (kw)  =  0  .  (12) 

u 

Equation  (12)  can  be  interpreted  either  as  the  source-free  Poisson  equation  in 
matter 


k^  e^(  cj)  <t>  (Ecj )  =  0 


(13) 


where  e  is  the  longitudinal  dielectric  constant,  or  as  the  solution  of  the  second 
order  differential  equation 


{j>  (xt)  +  « 


2 

P 


<t>  (xt)  =  0  , 


i.  e.  ,  the  oscillator  equation  for  the  hydrodynamic  oscillations  in  the  plasma. 
Returnirg  now  to  the  non-linear  equations,  we  first  solve  Eq.  (6)  for  i|j(lc<«))  and 
substitute  the  value  of  V(kco)  from  Eq.  (7)  to  obtain 


4<{kw) 


^  fa*  1 

— ^  <|>(kw)  /k  •  V(k-k',co-u)')(k-k') 

mw  w 


V(k'co' 


,  dk'dw' 


—  /k 

CO 


V(k'co')  ^  (k-k',a>-(o') 

(27r)^ 


(14) 


Further  substitution  of  Eq.  (14)  in  the  Poisson  equation  finally  gives 


5 


{1- 


CO  _  _  _ 

—2—  }  <j)  (koj)  =  — 2 — ~  /  ^  ■  V  (k-k',  w-M'){k-k') 


u 


CD 


(2<or 


4'iren  _  _  _ 

+  ^  °  j  Ic  •  V(k'w')  ij;  (k-k',  w-cd') 


dk'dcD' 

- ?■ 

(2cDr 


(15) 


To  bring  Eq.  (15)  in  a  more  lucid  form,  we  substitute  the  Vs  inside  the  integrals 
with  their  values  given  by  Eqs.  (6)  and  (7).  Thus 


CO 


k^  ct,(kco)  = 

OJ 


47r  en 


I 


k-  (k-k')  k-  (k-k')  , 


CO 


CO-CO 


(O’ 


<j>  (k-k',  co-co'  )(t)(k'co'  ) 


T-.  ,,dk'd< 


(Ztt) 


47r  e.  n  „  t; 


mco 


o  J  .j,  (■lk'(o')iJ; (k-k',  co-co') 


dk'dco' 

liT? 


(16) 


+  .  .  . 


The  terms  that  have  been  neglected  could  be  further  reduced  by  successive 
substitutions  to  give  higher  order  terms;  however,  the  algebra  is  prohibitive. 
In  any  case  we  can  clearly  see  that  Eq.  (16)  is  an  integral  equation  for  the 
potential  <j)  (kco)  and  its  right-hand  side  could  be  interpreted  as  a  source  term 
for  the  Poisson  equation  in  matter,  i.  e.  , 


6 


(17) 


u 


k  {1 - <|>  (kw)  =  47rp(koj) 


w 


Our  problem  now  is  to  find  a  solution  for  the  integral  equation.  Since  p(kw)  is 
expected  to  be  a  small  quantity,  it  becomes  obvious  that  we  should  seek  solu¬ 
tions  close  to  the  ones  given  by  the  linear  hydrodynamic  model.  According  to 
the  model,  density  oscillations  in  the  plasma  can  occur  at  only  one  frequency 
o)  and  at  an  arbitrary  wave  number  k.  We  know,  however,  from  the  solution 
of  the  B-V  equation  that  this  is  only  approximately  true;  i.e.  ,  the  oscillation 

frequency  is  given  by  w  =  oj  ,  k)  except  in  the  limit  of  k  ->■  0  where  w  =  «  . 

P  P 

Selecting  therefore  an  arbitrary  but  small  wave  number  k^,  we  assume  that  at 
a  first  approximation  the  potential  in  the  hydrodynamic  case  is  given  by 


<t)(kw)  =  <t)^  (Ztt  )  S(k-kp)S(a)-Wp)  , 


(18) 


where  <t>^  is  the  amplitude  of  the  potential.  However,  the  integral  equation  16 
expresses  the  fact  that  if  non-linear  terms  are  accovinted  for  there  are  contribu¬ 
tions  to  the  potential  from  all  k'  and  w'.  Stated  another  way,  non-linearities 
produce  a  slipping  of  the  frequency  and  wave  number,  a  process  very  similar 
to  incoherent  scattering.  In  fact,  if  a  plasma  oscillation  is  to  be  viewed  as  a 
density  wave  propagating  in  the  plasma,  then  it  becomes  clear  that  it  can  be 
scattered  out  of  the  direction  it  is  propagating.  As  the  wave  is  scattered  out 
of  the  forward  direction  it  loses  a  part  of  its  energy.  As  a  first  approximation 
we  assume  the  attenuation  to  be  small  and  separate  (}>(kw)  in  two  parts:  the 

coherent  part  consisting  of  the  propagating  density  wave  at  co  -  k  ,  and  the 

P  P 

incoherent  part  i|j^^^(kw): 


<j>{kw)  =  <t>^(27r  g(k-kp)  5(w-Wp)  +  <t)^^^(kw)  . 


(19) 
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We  introduce  now  Eq.  (19)  in  Eq.  (16)  and  obtain 


(1 - ^ — )  (l)^(27r)*^^(aj-a)^)  S(Ic-Ic^)  +  k  (1  — ^  )  (j)^  Nkw) 

u 


t-4  .. 

— ^  ())  f  k-k'  (Zjr)"^  S(w'-a)  )  S{k'-k  )  i|j(k-k',  co-oj') 

o  CO  P  P  (Ztt) 


,  ^p  r  cok'k'  j.(l)/T t  i\./T  Ti  i\  dk'dw'  , 

+  I  — i —  ‘I’  '(k'<o')i)j(k-k',  w-ca')  - +U  , 

CO  “  (2ir) 


(20) 


where  (T consists  of  terms  contributed  by  the  first  integral  in  Eq.  (16),  We  can 

see  by  inspection  that  these  terms  are  both  small  (since  they  are  proportional 

to  and  <t)^^^(|)^^^)  and  incoherent  with  respect  to  <j>^  and  (j>^^^(kw).  We 

will  therefore  neglect  them.  From  Eq.  (20)  we  clearly  have 

2  2 

(1  .  !^)4,(l)(kco)  =  <t>  - —  ^  fR-T:  ,  w-co  ) 


CO 


(0 


CO  p 

p 


CO  ^n(Ic-'k  ,  w-w  ) 

=  A  ^  -k-Tt  - E - E. 

^o  (0  p  n 


(21) 


and  therefore 


<b  coco  (k'K  )  On(H-k  ,  w-w  ) 

P  .  P - 2 - El 


,272  2, 

k  (<o  -cOp) 


(22) 


2  2 

where  the  division  by  co  -00^  will  be  interpreted  as  follows.  We  first  expand 


1 

T~T 

CO  -co_ 


^  (-J _ L_ 

2co„  w-co„  co+co„ 

p  p  p 


and  then  write 


8 


lim  - T— : — 

cr^  0  w-Wp  ^ 


1 

w-co 

p 


iri  £  (w-w  ) 
P 


(23) 


along  with  a  similar  expression  for  (co+cOp)  .  This  corresponds  to  the 
assumption  of  a  negligibly  small  loss  term  in  either  the  dielectric  constant 
e(w)  or  the  equation  of  motion  for  ^(xt).  Inserting  Eq.  (22)  back  in  Eq.  (20), 
we  obtain 

k^e(u)ct>^(27r  )^^(to-a)p)  £(k-kp)  = 

(k-k')(k-kp)  ^nCk'-iCp,  w'-ojp)  rfn(k-k',  w-w') 

772  ,2  t  2  ,,  .4 

k'  o)'  -o)  n  ( 27r  ) 

P  o 


This  equation  can  be  interpreted  as  an  equation  describing  the  scattering  of  a 

coherent  density  wave  at  w  ,  k  from  density  waves  traveling  at  directions  k' 

P  P 

and  frequencies  w'.  It  is  clear  that  the  density  deviations  rfnfk'-k  ,  w-w  ), 

P  P 

^  n(k-k',  w-w' )  inside  the  integral  of  Eq.  (24)  represent  density  waves  with 

respective  wave  numbers  and  frequencies  k'-k  ,  w-w  ;  k-k',  w-co'.  However, 

P  P 

it  is  also  clear  from  the  left-hand  side  of  Eq.  (24)  that  we  are  interested  in 


the  case  where  k  =  k  ,  w  =  w  ,  i.  e.  ,  in  the  forward  direction.  (See  Fig.  1.  ) 

P  P 


Fig.  1 


Density  Waves 


(24) 
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This  corresponds  to  imposing  the  conditions* 


(k'  -  kp)  +  (k  -  k')  =  0  i.  e.  ,  k  =  kp 


(25) 


(oj'  -  Wp)  +  (w  -  0)')  =  0 


1.  e.  ,  w  =  u> 


(26) 


(5) 

Taking  now  the  average  on  both  sides  of  Eq.  (24)  and  noting  that 


in(k'-k  ,  w'-w  )  gn(k-k',  cj-oj')  =(27r )  <  |  4n(k' -k  ,  w' -co  )  j  >5(k-k  )^(c»)-w  ) 
P  P  P  P  P  P 


and  that  <4>q>  =  <t*Q»  we  derive  from  Eq.  (24) 


k  [  e  (w)  +S  €(kw)]  (2ir  )  (}>  6(k-k  )  ( (w-co  )  =  0 

o.  p  p 


where 

^e(kcj)  = 


1  0)  (k  •  k')^ 

1-2  w  J  ,_,2 


CO 


,2  2 

CO'  -CO 


<]Mk'-k^, 
T 


P'  P'  ‘  dk'dco  ^27) 


n 


(27r) 


2 

where  now  <l^n|  >  is  the  mean  square  fluctuation  of  the  density.  We  note 
that  we  cannot  derive  information  concerning  the  amplitude  of  plasma  oscilla¬ 
tions  from  the  hydrodynamic  or  the  B-V  equations.  Nevertheless,  it  is  the 
thermal  motion  in  the  microscopic  scale  that  initiates  the  plasma  oscillations 
in  the  macroscopic  scale  and  therefore  the  amplitude  of  the  oscillations  can 
only  be  derived  from  thermodynamic  reasoning.  Appendix  A  is  devoted  to  the 
evaluation  of  the  fluctuation  spectrum  of  a  hydrodynamic  plasma  by  making  use 
of  Nyquist's  theorem.  We  find  that 


*  If  instead  of  density  waves  we  were  to  consider  particles,  Eqs.  (25)  and  (26) 
would  represent  momentum  and  energy  conservation. 
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(28) 


1  (Jn(kco)  1  ^ 


>  =  {S(w-Wp)  -  S(w+c^)} 


where  is  the  Debye  length.  Since  the  fluctuation  spectrum  has  just  two 

lines  at  ±  o)  and  no  width,  we  are  entitled  to  call  the  fluctuation  a  plasma 
P 

oscillation.  The  fluctuation  spectrum  derived  from  the  B-V  equations  has  a 
finite  width  which,  in  effect,  corresponds  to  Landau  damping  ;  however,  we 
can  readily  show  that  in  the  limit  of  0,  the  B-V  spectrum  is  again  given 

by  Eq.  (28).  With  the  available  information  we  can  now  compute  the  correction 
to  the  dielectric  coefficient  e(^oj).  Calling 

Se  =  Se'  +  iSe" 


K  =  k'  -  k  ,  -A  =  <*)'-  w 


and  using  Eq.  (23),  we  obtain 


2  2 

-  1  J!£_  1  _ ^ 


(29) 


L"  =  +4  ^  )- W))  ^  (30) 


n 


where 


Explicitly  we  have  for 

-  r  <"' '  ( 

J  —2  ^'2  ,0  ^3  ZFh~  3  0) 


D  2 


k  k' 


(27r  )■ 


(27r  )■ 


<jSn(K.n)!^>  =  ((O' -Si 


o 

where  g(k)  is  the  value  of  the  integral  between  0  and  some  k 

11 


max 


(32) 


The  dispersion  equation  then  becomes 


e  +  6€'  =  0  =  1  .  ^  +  g(k) 

CO  o 


which  gives 


;,2 

2  2  /-r. 

W  =  -  SW 


(0 


Stt  n  CO 
o 


2  /r\  D  2 

CO  -  g(k)  -= -  CO 

p  S'  '  37rn  p 


=  -^-g(k)(^  3^)  .  (33) 

^  o 

This  implies  that  the  frequency  decreases  because  of  non-linearities.  As  a 
result  the  group  velocity 

9(0 

'"g  "  ?K 

becomes  negative.  If,  however,  the  vector  k'  in  Fig.  1  were  to  be  reversed, 
the  sign  of  g(k)  in  Eq.  (33)  would  also  reverse  and  in  this  case 

2  2  ,  /T\  1  \ 

“  =  “p  +  g(k)  (_  3^) 

^  o 

As  a  result  v  >0. 

g 

From  Eq.  (30)  we  find  that  66"  =  0.  As  a  consequence  and  to  the 
approximation  we  are  using,  there  is  no  attenuation  in  the  propagating  density 
wave.  In  conclusion  we  may  say  that,  to  first  order,  non-linearities  affect 
only  the  dispersion  equation. 
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III. 


Non-Linear  Vlasov  Equation 


As  mentioned  in  the  previous  section,  the  hydrodynamic  model  of 
the  plasma  is  inadequate  in  the  sense  that  it  does  not  give  a  dispersion  re¬ 
lation  and  a  damping  term  in  the  linear  approximation.  It  becomes  therefore 
necessary  to  study  the  non-linear  Vlasov  equation  in  order  to  obtain  a  cor¬ 
rection  to  the  dispersion  relation  and  to  the  damping.  We  outline  here  the 
method  that  could  be  useful  in  a  program  of  this  nature.  However,  we  do 
not  attempt  to  give  an  answer  to  the  problem  because  the  equations  turn  out 
to  be  very  involved.  Instead  we  use  some  of  the  results  to  improve  the 
corresponding  hydrodynamic  formulation.  Vlasov's  equation  for  the  dis¬ 
tribution  function  f(x^)  has  the  form 

+  v'V-f+  —  Vd>-Vf=0  (34) 

ot  X  m  ^  V 

while  the  Poisson  equation  for  a  neutralized  electron  plasma  is 

V^<))  =  4jre  ^  fj^  (x^)  dv  (f  =  f^  +  fj^)  (35) 

Fourier  transformation  in  space  and  time  of  both  equations  gives 


(-icj  +  iR  ■  v)  f,  (Rwv)  +  i  —  R  •  V-  f  <j>(kw) 
'  '  1'  '  m  V  o^'  ' 


+i^  rc|,(k'a,')  R'-V-  f  (R-R'.u-a,')  =0 

^  ^  ^  (Zwr 


k  <j)(Rcj)  =  -  47re 
Solving  for  fj^(Rcov)i  we  obtain 


J  fj(Rwv) 


dv 


fj^(Rwv)  = 


k-  V-f 
e  V 


m 


V— I  _ 

— <j)(R<iO)  +  ^  \<))(k'w') 

co-R-  V 


CO  -  R  •  V 


(2^) 


(36) 
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We  integrate  Eq.  (36)  over  v  and  insert  in  the  Poisson  equation 


47re^  1 

r*  k.  V-f 

\  J  ^ 

mk^  " 

)  co-1c*7 

47re^  f 
m  J 

1  (t)(k'co') 

k' .  (k-k',  to-co',  v) 


v^r - dk’do;' 

co-k,  V  (2  7r)^ 


dv  . 


The  quantity  in  the  square  brackets  is  by  definition  the  longitudinal  dielectric 
coefficient  (kw),  and  therefore,  we  can  now  write  the  equivalent  of  Eq,  (16) 


k^£^(koj)<|)(koj)“ 


m 


k',V— f,  (k-k',  w-o)',  v)  , 

•cj') - -  - - -  dv 

u-k‘v  (2ff ) 


(37) 


=  4?:  p  (kcj) 

Equation  (37)  is  an  integral  equation  for  4>(kcu),  and  it  represents  the  same 
physical  problem  as  the  equivalent  hydrodynamic  equation.  It  presents, 
however,  a  mathematical  problem  of  larger  complexity  because  of  the 
denominator  (w-k' v).  To  bring  Eq.  (37)  in  better  correspondence  to  Eq.  (16), 
we  write  _  _  _  _ 


I 


k'*  Vvfi(k,A.v) 


dv  =  y 


k-  k'fj^(k,JX.,  v) 
(w-k*  v)^ 


dv 


w-k.  V 

It  then  becomes  obvious  that  for  k  -*■  0  we  obtain  Eq.  (16).  To  proceed  from 
Eq.  (37)  we  could  use  an  iteration  procedure,  or  we  could  write  as  in  Section  2 

(t>(ku)  =  5(k-kp)6  (w-w^)  +  (|)^^^(kw)  (38) 

2  2  2  2  2  2 
with  w  =  CO  +  (3  /c  T/m)  k  =  co  +  v  k  . 

p  '  p-i  p  o 
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Then 

k^e^(kw)  <j>^(27r  ^(k-k  )  ^(w-u  )  +  k^  e^(k«)  (j)^^^(kw) 


.2  k'  •  V-f ,  (k-k  ,  w-oj  ,  v) 

/  _E - ^  ^  p — SL-L  d? 


m 


cj-k  •  V 


tT  .  “  /  o_ 


w-k  ‘  V 


or  finally  (see  Sec.  II) 


k^e^(k(j)  <j)  {2Tr)^S{k-k  )  5'(w-a)  )  = 
o  P  P 


dco’ 


(27r) 


dv  , 


2-  [k  -V— f,(k'-k  ,v')][k'-  V— f,  (k-k',  w-w' ,  v)]  ,,,,  , 

(ii^)2  4,  f  -..  p„  1. 1 — p, . p  -  - :Wdv'  . 

^  °  k'  €  (k'co')  (w'-k'  •  v')  (w-k-  v)  {Ztt) 


(39) 


Our  intention  was  to  derive  an  equation  of  the  form 


k^[€^w)  +  ^€(k«)]  <t>  (27r  )^  S(k-k  )  g  (w-w  )  , 

Kj  P  Jr 


and  from  this  the  corrected  dispersion  relation  e'  (kw)  +  Se*  (kw)  =  0.  How¬ 
ever,  we  have  not  been  able  to  make  progress  along  these  lines.  As  a 
hybrid  solution  to  the  problem  we  could  use  the  left-hand  side  of  Eq.  (39) 
and  the  hydrodynamic  expression  for  the  right-hand  side.  From  Eq.  (27) 
we  would  then  have 


1  “p  1  <,l  n(k'-k  ^ 

s"(k“)  =  -  -mf-  i  - uftt:  - f - ^ —  ,T-x ' 


e^(k'w') 


n 


(27r) 
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or  since  from  the  fluctuation  theorem  we  have 


<  1  din(kh))  1  ^  kT  h. 

n_  TT  2  . _ ^ 


(lie  47rn 


D 


^11 


o 

L,,- 


1 

e^(k  w) 


) 


(kw) 


TTU 


e^{kw)  1  ^ 


Equation  (40)  becomes 


(i) 


d€(kb>)  =  - 


2  oi 


5 


k  •  k'  ‘ 

p 

771^ 


2  2 

e'^(k'ca')-ie"^(k'a)')  ^  e'‘^(KXV) 

|e^(k'  .  co-)l^  |e^(KA)|^ 


dk'dcj 

(2;r)^ 


(41) 


If  we  use  the  hydrodynamic  formula  for  the  fluctuations,  the  integration  over 
w'  can  be  performed  with  ease.  The  results  are  interesting  for  the  imaginary 
part  of  6e(kw): 


6e"(kw) 


I?"  J  k'‘=  |£V«')|‘'  "o  P  “  “-"p  P 


=  + 


‘°{k_  k')^(k'-kp)2 
_2-  "2 

k^k''^  (2 


e"^{k'  ,  2o>  ) 


dE-  7  ’  ^p'  ^  6"^(k',  0)  }  j;?  'ID  1 

tt)^  l€^(k',2w)l^  |e^(k',  0)1^  " 


(42a) 


(k  •  k')^(k'-k  )^  €"^(k',  2w  ) 

_  E  P  P'  r  P' 

-  L  rzTT^  [ 


k'<k 


k  k 


D 


|€-^(k',2wp)|‘-  |e"^(k',0) 


6"^(k',  0)  ,  1 

L/C-  OA I  2  ^  ^ 


(ii=(ii 


{42b) 


where  the  integration  or  summation  extends  from  k=0  to  k^^,  the  Debye  wave 
number.  We  note  now  that  the  fluctuations  of  the  electric  field  in  the  plasma 
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are  given  by 


<  1 E  (kw)  I  \ 


8  tcT 


(43) 


and  therefore,  the  energy  is 


C  (ko,)  =  =  si 

^  on  TTw  I  L„  vi 


e  (kw)  I 


Thus 


6e"  (koj) 


(k  •k')^(k'-k 


1  1 


k'<k 


k  k' 


D 


m  n 
P  o 


(44) 


where  the  summation  is  for  wave  numbers  smaller  than  the  Debye  wave  number. 
The  imaginary  part  of  6e  gives  a  measure  of  the  attenuation  of  the  coherent 
density  wave  because  of  non-linearities  in  the  system.  Equation  (44)  can  now 
be  compared  with  a  similar  equation  quoted  by  Pines.  ^  In  the  Pines  picture 
this  attenuation  is  expressed  in  terms  of  the  rate  of  change  of  plasmons  due 
to  non-linear  interactions.  In  addition  to  66"  we  could  also  compute  6e'  to 
derive  a  correction  for  the  B-V  dispersion  equation. 

In  conclusion,  we  may  say  that  in  the  hydrodynamic  model  of  a  plasma 
only  the  dispersion  equation  is  affected  by  non-linearities.  The  indication, 
however,  is  that  in  the  B-V  model  of  a  plasma  non-linearities  introduce  not 
only  a  correction  to  the  dispersion  equation  but  also  an  attenuation. 
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Appendix  A;  Calculation  of  the  Density  Fluctuation 


If  an  external  potential  <j)  is  included  in  the  linearized  version  of 
the  hydrodynamic  equations,  we  have  from  Eqs.  (9-11) 


-icoi|j(R<o)  +  iR  *  V  (Rco)  =  0 

(A.l) 

“ico'V  (Rco)  +  iR  ^  <j)(Rco)  +  iR  ^  (t)®^^(Rco)  =  0 

(A.  2) 

-k^<j)(Rco)  +  47r  en^ijj{Rco)  =  0 

(A.  3) 

Multiplying  (A..  2),ty  R,  introducing  the  result  in  (A.  1)  and  making  use  of  the 
Poisson  Eqnatipn  we  have  successively 


4i(koj)  = 


:  •  V  (Rw) 

- 2 - 

CO 


1,2  ,  2  ^ 

e  K  e  k  ,  ext.r  . 

~m  ""2  m 


=  ilS^  + 


mco 


m 


CO 


Thus 

(co^-«p)  i1j(R«)  =  —  k^(l)®^^Rw) 
or 

CO  -  CO 

p 

i|j(kco)  is  the  linear  response  of  the  hydrodynamic  plasma  to  an  external  poten¬ 
tial  Since  iIj(Rco)  =  6n(kco)/n^,  we  have  further 

2^2  2 
p(Rco)  =  e6n(kw)  =  ^ 

CO  -  CO 

p 
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2  2 

where  e  =  1-  w  /w  .  The  current  is  given  by 
P 


j(kw)  = 


.  CO 

ICO  p  —  ext  .7  V 

-r--z  (>“1 

U>  -  O) 

p 


(A.  5) 


^"©xt "  ““  ext 

where  E  (ku)  =  -  ik<t>  (kw).  From  Nyquist's  theorem  we  know  that  the 


mean  square  fluctuation  for  j  is 


CJ 


W  -  W 

P 


(A.  6) 


and  therefore 


<  j6n(koj)  1^  >  =  ^-*2  <  j  j(kw)  > 

ca)  e 


If  we  interpret  the  singularity  of  (A.  6)  as 
2 


^  {6(«.«  )  -6(w+«  )}. 

CO  -  U)_  ^  P  P  tr  f 


we  readily  obtain 


6n(kcj) 


kT 


n 


47re^n 


k^co 

.{6(w-Wp)  -6(w+Wp)  } 


(A.  7) 


which  integrates  to 


J'  <  j  5n(kco)  \  ^>  dco  = 


.^4 


where 


Xq  =  KT/47re^n^ 


is  the  Debye  length. 
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